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ABSTRACT 

Cygnus X-l is a binary star system that is comprised of a black hole and a 
massive giant companion star in a tight orbit. Building on our accurate distance 
measurement reported in the preceding paper, we first determine the radius of the 
companion star, thereby constraining the scale of the binary system. To obtain 
a full dynamical model of the binary, we use an extensive collection of optical 
photometric and spectroscopic data taken from the literature. By using all of the 
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available observational constraints, we show that the orbit is slightly eccentric 
(both the radial velocity and photometric data independently confirm this result) 
and that the companion star rotates roughly 1.4 times its pseudosynchronous 
value. We find a black hole mass of M = 14.8 ± 1.0 M , a companion mass of 
M opt = 19.2 ± 1.9 M Q , and the angle of inclination of the orbital plane to our line 
of sight of % = 27.1 ±0.8 deg. 



Introduction 



During the past 39 years, many varied estimates have been made of the mass M of 
the black ho le in Cygnus X-l. A t one extreme, acting as a devil's advocate against black 
hole models, iTrimble et al.l ( 119731 ) proposed a model based on a distance D ~ 1 kpc which 
gave a low mass of M < 1 M , suggestive of a neutron star or white dwarf, not a black 
hole. Several other low-mass models are summarized, considered, and found wanting by 
Boltonl ( 119751 ). By contrast, all conventional binary models that assume the secondary 
companion is a massive O-type supergiant find a large — but uncertain — mass for the compact 
object that significantly exceeds the maximum stable mass for a neutron star of ~ 3 M G 
( Kalogera fc Bavm 19961), he nce requiring a black hole. For example, using geometrical 
arguments, iPaczyhskil ( 119741 ) computed the minimum mass for the compact object as a 
function of the distance and fo und M > 3.6 M Q for D > 1.4 kpc. Based on dynamical 
modeling iGies fc Boltonl (I1986T) found M > 7 and a probable mass of M opt = 16 M Q 
for the companion star, and iNinkov et al.l (119871 ) found M = 10 ± 1 M, 



q (by assuming 



M, 



opt 



20 M, 



e) 



However, these mass estimates, and all such estimates that have been made to date, 
are very unc ertain because th ey are based on unsatisfactory estimates of the distance to 
Cygnus X-l ( iReid et al.l 120111 ). The str ong effect of distance on the model parameters is 
obvious from an inspection of Table 4 in ICaballero-Nieves et al.l (120091) (and also Tabl e 1 of 
Paczyhskil 1 19741 ) . For the dynamical model favored by ICaballero-Nieves et al.l ( 120091 ). and 
over the wide range of distances they consider, 1.1-2.5 kpc, the radius of the companion 
star and the mass of the black hole are seen to vary by factors of 2.3 and 11.7, respectively. 
Thus, in order to obtain useful constraints on the system parameters, it is essential to have 
an accurate value of the source distance, as we have demon strated for two ext ragalactic black 
hole systems that contain O-type supergiants, M33 X-7 (jOrosz et al.l 120071 ) and LMC X-l 
( jOrosz et al.l 120091 ). whose distances are known to several percent accuracy via the cosmic 
distance ladder. 



In this paper, we use a distance from a trigonometric parallax measurement for Cygnus 
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X-l (IReid et al.ll201lf ). which is accurate to ±6%, and previously-published optical data to 
build a complete dynamical model of the Cygnus X-l binary system. Not only are we able 
to strongly constrain the principal parameters of the system, we are also able to obtain the 
first constraints on the orbital eccentricity, e, and the deviation of the period of rotation of 
the companio n star from the or bital period, characterized by the non-synchronous rotation 
parameter, Q ( lOrosz et al.ll2009l ). Of principal interest are our precise determinations of the 
black hole mass M and the orbital inclination angle i. As we show in the paper that follows 
( iGou et al.ll201ll ). it is our accurate values for the three parameters M, i and D that are the 
key to determining the spin of the black hole. 



2. Dynamical Modeling 

The mass of the black hole can be easily determined once we know both its distance 
from the center of mass of the O-star and the orbital velocity of the star. Since optical 
spectroscopy only gives us the radial component of velocity, we must also determine the 
inclination of the orbital plane relative to our line of sight in order to infer the orbital 
velocity. Furthermore, since the star orbits the center of mass of the system, we must also 
obtain the separation between the two c omponents. We determin e the needed quantities 



using our comprehensive modeling code (lOrosz fc Hauschildtl |2000| ). Our model, which is 



underpinned by our new measurement of the distance, makes use of all relevant observational 
constraints in a self-consistent manner. We discuss details of the data and modeling below. 



2.1. Optical Data Selection 



There is no shortag e of published observational data for Cygnus X-l in the literature. 
Brocksopp et al.l (Il999al ) provide U, B, and V light curves containing nearly 27 years worth 
of observations (1971-1997) from the Crimean Laboratory of the Sternberg Astronomical 
Institute. The binned light curves contain 20 points each and are phased on the following 
ephemeris: 

Min I = JD 2441163.529(±0.009) + 5.599829(±0.000016)£, (1) 

where Min I is the time of the inferior conjunction of the O-star, and E is the cycle number. 
In addition, iBrocksopp et al.l (ll999al ) provide 421 radial velocity measurements, which are 
phased to the same ephemeris. The light and velocity curves were kindly sent to us by C. 
Brocksopp. 



In addition to the velocity data of Broscksopp et al., we also made use of the radial 
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velocities published by IGies et al.l (120031 ) . We fitted a sine curve to their data and those of 
Brocksopp et al. (1999a). Fixing the uncertainties on the individual velocity measurements to 
7.47 km s _1 and 5.0 6 km s -1 , respecti vely (these values give \ 2 ~ N), we found K = 74.46 ± 
0.51 km s _1 for the iBrocksopp et al.l data and K = 75.57 ± 0.70 km s _1 for the IGies et al. 



data. The la intervals of the respective i^-velocities overlap, the residuals of the fits show 
no obvious structure, and there is no notable difference between the residuals of the two 



data sets, apart from the slightly greater scatter in the IBrocksopp et al.l data. We therefore 
combined the two sets of radial velocities while re moving the respective systemic velocities 
from the individual sine curves. We phased the IGies et al.l data on the above ephemeris 
before merging these data with the velocity data of Broscksopp et al. The combined data 
set has 529 points. 

All of these velocity and light-curve data are discussed further and analyzed in Sec- 
tion 2.4, where they are shown folded on the orbital period and seen to exhibit minimal 
scatter about the model fits. This small scatter in these data sets, each spanning a few 
decades, attests to the strongly dominant orbital component of variability. Meanwhile, 
Cygnus X-l is well-known to be variable in the radio and X-ray bands, including major 
transitions between hard and soft X-ray states (see Figure 1 in Gou et al. 2011). This raises 
the question of whether any non-orbital variability in the light-curve data could significantly 
affect the component masses and other parameters determined by our model. (We focus on 
the light-curve data, which are more susceptible to being affected by variability.) 

We believe that our results ar e robust to such non-orbi tal variability for several reasons, 
including the following: (1) While IBrocksopp et al.l (Il999bl ) note that the U, B, and V light 
curves were correlated with each other, they found no correlation with the radio and X-ray 
light curves, which is reasonable given that the time-averaged bolometric X-ray luminosity 
is only < 0.3% of the bolometric luminosity of the O-star (Section 2.3). (2) Although 
Brocksopp et al.l (ll999aT) did not specifical ly discuss the photometric variability seen in their 
27 year data set, IBrocksopp et al.l (Il999bl ) do discuss the multiwavelength variability of the 
source over a 2.5 year period. They noted that there is very little variability in the optical 
light curves apart from the dominant ellipsoidal/orbital modulations (see Section 2.4) and 
one weaker modulation with about half the amplitude on a 142 day period (which is thought 
to be related to the precession of the accretion disk). (3) The light-curve and velocity data 
give remarkably consistent results for the small (but statistically very significant) measured 
values of eccentricity and argument of periastron (Section 2.4.2). (4) As highlighted above, 
the small scatter in the folded light curves for a data set spanning 27 years is strong evidence 
against a significant component of variability on time scales other than the orbital period. 
In summary, we conclude that non-orbital variability is unlikely to significantly affect our 
results. 
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2.2. Stellar Radius, Temperature, and Rotational Velocity 



In order to constrain the dynamical model, it is crucial to have a good estimate of 
the radius of the companion star. However, customary methods of determining this ra- 
dius fail because the Cygnus X-l system does not exhibit eclipses nor does the compan- 
ion star fill its Roche equipotential lobe. We obtain the req uired estimate of the stellar 
radius as we have done previously in our study of LMC X-l (jOrosz et al.ll2009l ). The ra- 
dius, which critically depends on distance, additionally depends on the apparent magni- 
tude of the O-type star and interstellar extinction, and also on the effective stellar tem- 
perature and corresponding bolometric correction. The absolute magnitude of the star is 
Mibs = K+BCxiTeg, <?) — (5 log-D — 5)— O.llAy, where K is the apparent .fT-band magnitude, 
BCk is the bolometric correction for the if-band, D is the distance, and Ay is the extinction 
in the V"-band. The luminosity and radius of the star in solar units are L = io _0 - 4 ( Mabs ~ 4 - 71 ) 
and R = \/L(577Q/T ^f)^, respectively. In computing these quantities, we u se D = 1.86ln I? 
kpc (IReid et al.ll201ll ) and a if-band apparent magnitude of K = 6.50 ±0.02 (ISkrutskie et al. 



20061 ). which minimizes the effects of interstellar extinction. For the if-band extinction, we 
adop t E(B-V) = 1.11 ±0.03 and R v = 3. 02 ±0.03 (e.g.. Ay = 3.35 jCaballero-Nieves et al" 



20091 ) and use the standard extinction law (ICardelli et al.lll989l ). The bolometric corrections 
for the if -band were com puted using the OSTAR2002 grid of models with solar metallic- 
ity (ILanz fc Hubenyl 120031 ). We note that the if-band bolometric corrections for the solar 
metallicity models and the models for half-solar metallicity differ only by 0.02 dex (T. Lanz 
private communication), so our results are not sensitive to the metallicity. 

Figure [1] shows the derived radius and luminosity of the star as a function of its assumed 
temperature in the range 28,000 < T cff < 34,000 K. For T eS = 28,000, the radius is 
R dist = 19.26 ±0.98^(11 As the temperature increases, the radius decreases rapidly at first, 
and then it plateaus midway through the range, attaining a value of -Rdist = 16.34 ± 0.84 R Q 
at T c fj = 34, 000 K. Meanwhile, the luminosity increases with temperature, rising from 
L = 2.1 x 10 5 L G , at T Pff = 28, 000 K to L = 3.2 x 10 5 L G , at T cff = 34, 000 K. 



The effective temperature of t he companion star can be determined from a detailed anal- 

ysisof UV and optical line spectra ( Herrero et alJ[l995 ; Karitskaya et al.ll2005 ; Caballero-Nieves et al 
20091 ). However, it is often difficult to determine a precise temperature for O-type stars ow- 
ing to a correlation between the effective temperature T e g and the surface gravity parameter 
logg. Model atmospheres with slightly smaller values of T e g and logg give spectra that are 
very similar to those obtained for slightly higher values of these parameters. Fortunately, 



1 We use the notation -Rdist to denote the stellar radius derived from the parallax distance, and i? pt to 
denote the stellar radius derived from the dynamical model. 
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log g for the companion star in Cygnus X-l is tightly constrained for a wide range of assump- 
tions about the temp erature, mass r atio, and other parameters because of a peculiarity of 
Roche- lobe geometry (jEggletonlll983l ). We show below that our dynamical model constrains 
the value of the surface gravity to lie in the range log g = 3.30 — 3.45 (where g is in cm s~ 2 ). 



Based on an analysis of both optical and UV spectra, ICaballero-Nieves et al.l (120091 ) 
obtained for their favored model T e ff = 28,000 ± 2,500 K and a surface gravity log g > 
3.00 ± 0.25, which is outside the range of values implied by our model. Based on the plots 
and tables in ICaballero-Nieves et al.l (120091 ). we estimate a best-fitting temperature of T cS = 
30, 000 ± 2, 500 K when the surface gra yity is forced to l i e in t he range determined by our 
dynamical model. Using optical spectra, iKaritskaya et al.l ( 120051 ) found T e ff = 30, 400 ±500 K 
and log q = 3. 3 1 ±0.0 7, which is consistent with our dynamically-determined value. Likewise, 
Herrero et al.l ( 119951 ) report T e g = 32,000 and log g = 3.21 (no uncertainties are given in 



their Table 1). In the following, we adopt a temperature range of 30, 000 < T cS < 32, 000 K. 
After considering several previous determinations of the proje cted rotational velocity 



of the O-type star and corrections for macroturbulent broadening, ICaballero-Nieves et al. 
(120091 ) adopt V TOt sini = 95 ± 6 km s _1 . We use this value as a constraint on our dynamical 
model, which we now discuss. 



2.3. ELC Description and Model Parameters 



The ELC (Eclipsing Light Curve) model (lOrosz fc Hauschildtl |2000| ) has parameters 
related to the system geometry and parameters related to the radiative prop erties of the star. 
For th e Cygnus X-l models, the orbital period is fixed at P = 5.599829 days (IBrocksopp et al. 



1999a]). Once the values of P, the K- velocity of the O-star and its mass M opt are known, the 
scale-size of the binary (e.g., the semimajor axis a) and the mass of the black hole M are 
uniquely determined. With the scale of the binary set, the radius of the star R opt determines 
the Roche-lobe filling fraction p. Not all values of R opt are allowed (for a given P, M opt , and 
K): If R opt exceeds the effective radius of the O-star's Roche lobe, we then set the value of 
p to unity. 

The main parameters that control the radiative properties of the O-star are its effective 
temperature T eS , its gravity darkening exponent (3, and its bolometric albedo A. Following 
standard practice for a star with a radiative outer envelope, we set /3 = 0.25 and A = 1. 

The ELC model can also include optical light from a flar ed accretion disk. In the cas e 
of Cygnus X-l, the O-star dominates the optical and UV flux (ICaballero-Nieves et al.ll2009l ). 
where the ratio of stellar flux to accretion disk flux at 5000 A is about 10, 000. Consequently, 
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we do not include any optical light from an accretion disk. 

We turn to the question of the X-ray heating of the supergiant star and its effect on th e 



binary model. The X-ray heating is computed using the technique outlined in IWilsonl ( 119901 ). 
The X-ray source geometry is assumed to be a thin disk in the orbital plane with a radius 
vanishingly small compared to the semimajor axis (this structure should not be confused 
with much larger accretion disk that potentially could be a source of optical flux). Points on 
the stellar surface "see" the X-ray source at inclined angles, and the proper foreshortening 
is accounted for. 

Measurement of the broadband X-ray luminosity of Cygnus X-l (L x boi; hereafter in units 
of 10 37 erg s -1 , adjusted to the revised distance of 1.86 kpc) requires special instrumentation 
and considerations. There are soft and hard states of Cygnus X-l (e.g., Gou et al. 2011). 
The hard state is especially challenging because the effective temperature of the accretion 
disk is relatively low (T < 0.5 keV), while the hard power-law component (with photon index 
~ 1.7) must be integrated past the cutoff energy ~ 150 keV (Gierlinski et al. 1997; Cadolle 
Bel et al. 2006). Since the ground based observations (i.e., photometric data and radial 
velocity measurements) are distributed over many years, both the range and the long-term 
average of the X-ray luminosity must be estimated. 

The archive of the Rossi X-ray Timing Explorer (RXTE) contains several thousand 
observations of Cygnus X-l collected in numerous monitoring campaigns conducted over 
the life of the mission. We processed and analyzed 2343 exposure intervals (1996 January 
to 2011 February; mean exposure 2.2 ks) with the PCA instrument, and we then computed 
normalized light curves in four energy bands (2-18 keV) in the manner described by Remillard 
& McClintock (2006). In the hardness- intensity diagram, the soft and hard states of Cyg 
X-l can be separated by the value of hard color (HC; i.e., the ratio of the normalized PCA 
count rates at 8.6-18.0 versus 5.0-8.6 keV), using a simple discrimination line at HC = 0.7. 
On this basis, we determine that Cygnus X-l is found to be in the hard state 73% of the 
time. 

Zhang et al. (1997) studied the broadband spectra of Cygnus X-l with instruments of 
RXTE and the Compton Gamma Ray Observatory (CGRO). During the hard and soft states 
of 1996, they found broad-band X-ray luminosities in the range 1.6-2.2 for the hard state 
and 2.2-3.3 for the soft state. Samples of the hard and intermediate states during 2002-2004 
with the International Gamma-Ray Astrophysics Laboratory (INTEGRAL) yielded L x boi in 
the range 1.2-2.0 (Cadolle Bel et al. 2006). Additional measurements of the soft state with 
Beppo-SAX (Frontera et al. 2001) found L x boi m the range 1.7-2.1, while Gou et al. (2011) 
analyzed bright soft state observations by ASCA and RXTE (1996) or Chandra and RXTE 
(2010) to determine L x t, i in the range 3.2-4.0. 
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To get a rough idea how these special observations relate to typical conditions, we used 
the available contemporaneous RXTE observations to scale the measured L xho \ values against 
the PCA count rates, considering hard and soft states separately. We then estimate that the 
average hard and soft states would correspond to L x boi values of 1.9 and 3.3, respectively. 
Finally, the time-averaged luminosity would then be roughly 2.1 x 10 37 erg s -1 , which cor- 
responds to 0.01 of the Eddington limit for Cygnus X-l. This value is much smaller than 
the bolometric luminosity of the O-star (Lr \ fh 7.94 x 10 39 erg s _1 ), and so we expect that 
X-ray heating in Cygnus X-l will not be a significant source of systematic error. We ran 
some simple tests in which we increased the bolometric X-ray luninosity by up to an order 
of magnitude, and found the light curves to be essentially identical. 

When computing the light curve, ELC integrates the various intensities of the visible 
surface elements on the star. At a given phase and viewing angle, each surface element on the 
star has a temperature T, a gravity log g, and a viewing angle /J = cos 6, where 9 is the angle 
between the surface normal and the line-of-sight. ELC has a pre-computed table of specific 
intensities for a grid of models in the T e g — logg plane. Consequently, no parameterized 
limb darkening law is needed. For the specific case of Cygnus X-l, the range of temperature- 
gravity pairs on the star contained point s that are outside the OSTAR2002+BSTAR2006 
model grids ( jLanz fc Hubenyl 12003k 120071 ) . We therefore computed a new grid of models, 
assuming solar metallicity. The range in effective temperature is from 12,000 K to 40,000 
K in steps of 1000 K, and the range in the logarithm of the surface gravity logg (cm s -2 ) 
is from 2.00 to 3.75 in steps of 0.25 dex. The wavelength-dependent radiation fields 
as a function of the cosine of the emergent angle /i, were computed from 232 spherical, 
line-blanketed, LT E models using the gener alized model stellar atmosphere code PHOENIX, 
version 15.04.00E (lAufdenberg et al.lll998l ). These models take into consideration between 
1,176,932 and 2,296,243 atomic lines in the computation of the line opacity. The radiation 
fields were computed at 26,639 wavelengths between 0.1 nm and 900,000 nm and at 228 angle 
points at the outer boundary of each model structure. The models include 100 depth points, 
an outer pressure boundary of 10~ 4 dyn cm -2 , an optical depth range (in the continuum at 
500 nm) from 10 -10 to 10 2 , a microturbulence of 2.0 km s _1 , and they have a radius of 18 Rq. 
(At the gravities of interest, the models are very insensitive to the precise value of the stellar 
radius.) A table of filter-integrated specific intensities f or the Johnson U, B , and V filters 
for use in ELC was prepared in the manner described in lOrosz fc Hauschildtl ( 120001 ) . 



As noted above, the obse rvational data we model include the U, B, and V light c urves 
from IBrocksopp et al.l ( 1l999ai ). and the radial velocities from iBrocksopp et al.l (jl999al ) and 
Gies et al.l (120031 ) . which are combined into one set. We use the standard x 2 statistic to 
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measure the goodness-of-fit of the models: 

2 fUi — U mo d{a)\ 2 ^ / Bi — B mod (a)\ 2 ^ fVi — V mo d(a) x ~ 



+ ^fRV l ~RV mod (a) x 2 

i=i 



(2) 



Here, the notation C/j refers to the observed U magnitude at a given phase i; Oi is the 
uncertainty; and U mod (a) is the model magnitude at that same phase, given a set of fitting 
parameters a. Similar notation is used for the B and V bands, and for the radial velocities. 

Initially, we assumed a circular orbit and synchronous rotation, which left us with five 
free parameters: the orbital inclination angle i, the amplitude K of the O-star's radial 
velocity curve (the "i^-velocity"), the mass of the star M opt , the radius of the star R opt , and 
a small phase shift parameter (p to account for slight errors in the ephemerides. We quickly 
found, upon an examination of the post-fit residuals, that circular /synchronous models are 
inadequate and more complex models are required. In the end, we computed four classes of 
models: (i) Model A has a circular orbit and synchronous rotation, and it uses the five free 
parameters discussed above: a = (i, k, M opt , R opt , 0). (ii) Model B has a circular orbit and 
nonsynchronous rotation, requiring one additional free parameter Q, which is the ratio of the 
rotational frequency of the O-star to the orbital frequency: a = (i, k, M opt , R opt , 4>, fi). (iii) 
Model C has an eccentric orbit and synchronous rotation at periastron. For this model there 
are seven free parameters, namely the five parameters for Model A plus the eccentricity 
e and the argument of periastron u: a = (i,k,M opt ,R opt ,<p,e,cu). Finally, (iv) Model D 
incorporates the possibility of nonsynchronous rotation and an eccentric orbit and has a 
total of eight free parameters: a = (i, k, M opt , R op t, 0, e, u>, Q). 

We also have additional observed constraints that do not apply to a given orbital phase, 
but rather to the system as a whole. These include the stellar radius (computed from the 
parallax distance i?dist for a given temperature) and the rotational velocity of V rot sin i = 
95 ± 6 km s _1 . For Cygnus X-l, X-ray eclipses are not observed, and hence there is no 



eclipse constraint. As discussed in lOrosz et al.l (120021 ). the radius and velocity constraints 
are imposed by adding additional terms to the value of x 2 : 



2 l -Ropt - -Rdist \ 2 . / Vrot sin % - 95 x 2 



XL = { aRd ) + { 6 ) ■ ( 3 ) 
Finally, for a given model, we have as our ultimate measure of the goodness-of-fit 

Xtot Xdata Xcon" (4) 
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As with any fitting procedure, some care must be taken when considering relative weights 
of various data sets. In the case of Cygnus X-l, the quality of the light curves is similar in all 
three bands. We therefore scaled the uncertainties by small factors in order to get x 2 ~ N 
for each set separately. This resulted in mean uncertainties of 0.00341 mag for U, 0.00155 
mag for B, and 0.00226 for V. The uncertainties for the radial velocity measurements are 
given in Section 2.1. 

2.4. ELC Fitting 

2.4-1- Principal Results 

The fits to the light-curve and radial- velocity data for Model A (left panels) and Model 
D (right panels) are shown in Figure 2, and the best-fit values of the parameters are given for 
all four models in Table 1. These results are for T cff = 31,000 K; similar figures and tables 
for other temperatures are presented in the Appendix. For a given temperature, the tabular 
data reveal a key trend: The total x 2 of the fit rapidly decreases as one goes from Model A to 
Model D, indicating that eccentricity and nonsynchronous rotation are important elements, 
and we therefore adopt Model D. For final values of the fitted and derived parameters given 
in Table 1, we take the weighted average of the values derived for each of the temperatures 
in the range 30, 000 < T cff < 32, 000 K. A schematic diagram of the binary based on our 
best-fitting model is shown in Figure 3. 

2.4-2. Need for N on- Zero Eccentricity 

In our dynamical model, the orbit is not quite circular and the stellar rotation is not 
synchronous at periastron. Although eccentric orbits based on light curve models have been 
considered in the past (e.g. Hutchings 1978; Guinan et al. 1979), most of the recent work 
considers circular orbits and synchronous rotation (e.g. Caballero-Nieves et al. 2009). While 
the eccentricity we find is small (e = 0.018±0.002 for Model D at T cff = 31000 K), it is highly 
significant. Allowing the orbit to be eccentric results in an improvement to the x 2 values 
of all three light curves and the radial velocity curve as well. As a separate check on our 
results, we fitted the light curves and the velocity curve separately and found the best-fitting 
values of e and u for each. For the light curves, which are derived from a homogeneous data 
set containing more than 800 observations spanning 27 years, we find e = 0.0249 ± 0.003 
and u = 305 ± 7°. For the velocity curve, we find e = 0.028 ± 0.006 and u = 307 ± 13°. In 
addition to the remarkable consistency of the results, note that the argument of periastron 
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we find is not consistent with e ither 90° or 270°, which indicates that the eccentricity is not 



an artifact of tidal distortions ( Wilson fc Sofialll976r lEatonll2008l ) 



2.4-3. Sensitivity of the Results on the Model Assumptions 

A glance at Table 1 shows some of the parameters, such as the inclination, seem to 
change considerably among the four models. These differences can be explained in part by 
how much weight one places on the external constraints, namely the radius found from the 
parallax distance and the observed rotational velocity. If we assume synchronous rotation, 
then the observed radius of 16.42 ± 0.84 R Q (for T eS = 31, 000) and the projected rotational 
velocity of V rot sini = 96 ± 6 km s _1 give an inclination of i = 39.8 ± 3.9°. We ran a model 
C that is identical to model C except we fixed the errors on the external constraints to be 
a factor of 10 smaller (e.g. 16.42 ± 0.084 R & and V^ ot sin2 = 96 ± 0.6 km s _1 ). We found 
an inclination of i = 39.6 ± 0.4°, which is consistent with the expected value. However, the 
fits are much worse, with x 2 — 629.4 compared to x 2 = 614.8. Likewise, we ran a model A' 
that is identical to model A except for the use of the same hard external constraints, and we 
similarly found i = 39.6 ± 0.4° and x 2 = 673.2. Thus, in this limiting case where the radius 
and rotational velocity are forced to their observed values, the inclination is determined by 
these quantities and does not depend on the eccentricity. 

In principle, the light curves should be sensitive to nonsynchronous rotation since the 
potential includes a term containing a fac tor of Q 2 , where Q is the ratio of the stellar rota- 



tional frequency to the orbital frequency (jOrosz fc Hauschildtl |2000| ) . As shown in Table 1 



we obtain improved fits to the light curve when the rotation is allowed to be nonsynchronous: 
model B has a smaller x 2 than model A (although both of these models are rejected since 
the eccentricity is nonzero), and model D has a smaller x 2 than model C. In the case of 
model D, note that the values of x 2 f° r each of the U, B, and V light curves has decreased 
relative to the values for model C. When the rotation of the O-star is not synchronous, the 
mapping between the radius, the inclination, and the observed projected rotational velocity 
is of course altered. Thus in model D we are able to satisfy the constraints imposed both by 
the observed radius and the observed rotational velocity while providing optimal fits to the 
individual light curves. 

Finally, we note that it is possible to achieve optimal fits to the light curves for an ec- 
centric orbit and synchronous rotation if we ignore the constraints provided by the measured 
values of the radius and rotational velocity. However, without these constraints there are 
degeneracies among solutions. For example, two solutions with similar x 2 values exist: one 
solution has i = 21.4°, R 2 = 22.98 .Rq, V mt smi = 98.4 km s~\ and x 2 = 597.0, whereas 
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another solution has i = 25.4°, R2 = 16.2 Rq, V TO tsmi = 63.0 km s _1 , and y 2 = 596.5. In 
the former case, the derived rotational velocity is consistent with the observed value, but 
the derived stellar radius is much too large. In the latter case, the derived stellar radius is 
consistent with the radius computed from the parallax distance; however, the derived rota- 
tional velocity is much too small. As mentioned above, allowing the stellar rotation to be 
nonsynchronous allows us to satisfy all of the constraints while optimally fitting the light 
curves. 



3. Discussion and Summary 

We find masses of M opt = 19.16 ± 1.90 M and M = 14.81 ± 0.98 M for the O-star 
and black hole, respectively. These estimates are considerably more direct and robust than 
previous ones, owing largely to the new parallax distance. This secure trigonometric distance 
was used to derive a radius for the O-star using the observed i^-band magnitude to minimize 
uncertainties due to extinction and metallicity. The derived O-star radius strongly constrains 
the dynamical model. In addition, we have also used the observed rotational velocity of the 
O-star as an observational constraint, and sought models that simultaneously satisfy all of 
the constraints. We have explored possible sources of systematic error by considering a wide 
range of possible temperatures for the companion O-star and circular and eccentric models 
separately. The optical light curves we used come from a long-term program which used 
the same instrumentation and should represent the mean state of the system quite well. 
Overall, the uncertainties in the dynamical model are as about as small as they can be in a 
noneclipsing system. 



111 



By way of comparison with previous work, a recent summary of mass estimates is given 



Caballero-Nieves et al.l (120091 ) . The four estimates given for the black hole mass that are 
based on dynamical studies (see Section 1) are quite uncertain and generally consistent with 
our result. Other estimates, whic h are based on less certain methods, generally imply a lower 



mass, M w 10 M . For example, lAbubekerov. et al.l ( 120051 ) repo rt a mass in the range 8.2 < 



M < 12.8 M & from their analysis of optical spectral line profiles. IShaposhnikov fe Titarchuk 



( 120071 ) found a mass in the range 7.9 < M < 9.5 M using X-ray spectral/timing data and 
a scaling relation based on the dynamically-determined masses of three black holes. Such 
methods, while intriguing, are not well established and are prone to large systematic errors 
compared to the time-tested and straightforward approach that we have taken. 

Our improved dynamical model for Cyg X-l enables other studies of this key black- 
hole binary. Using our precise measurement of the distance (Reid et al. 2011), we are now 
able to compute the stellar luminosity as a function of the assumed temperature, which 
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allows one to place the star on a H-R diagram with some confidence. We have also provided 
precise values of the component masses, the eccentricity, and the degree of nonsynchronous 
rotation, quantities which may be used to test binary evolutionary models (such a study will 
be presented elsewhere). Finally, with our accurate values for the distance, the black hole 
mass, and the orbital inclination angle, one can model X-ray spectra of the source in order 
to measure the spin of the black hole primary. In our third paper in this series, Gou et al. 
(2011), we show that Cyg X-l is a near-extreme Kerr black hole. 

We thank Catherine Brocksopp for sending us the optical light curves. The work of 
JEM was supported in part by NASA grant NNX11AD08G. This research has made use of 
NASA's Astrophysics Data System. 
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Table 1. Cygnus X-l Model Parameters 



parameter 


Model A 


Model B 


Model C 


Model D 


adopted 


i (deg) 


51.45 ± 7.46 


67.74 ±2.83 


28.50 ±2.21 


27.03 ± 0.41 


27.06 ± 0.76 


n 


1.000 


0.674 ± 0.043 


1.000 


1.404 ± 0.099 


1.400 ±0.084 


e 






0.025 ±0.003 


0.018 ±0.002 


0.018 ± 0.003 


u (dcg) 






303.5 ±5.1 


308.1 ± 5.5 


307.6 ±5.3 


Mopt (M©) 


20.53 ± 2.00 


26.27 ±2.41 


25.17 ±2.54 


18.97 ±2.15 


19.16 ± 1.90 


M (M©) 


7.37 ± 1.19 


6.98 ± 0.39 


15.83 ± 1.80 


14.69 ± 0.70 


14.81 ±0.98 


-Ropt (-Rq) 


15.07 ± 1.26 


16.41 ±0.72 


18.46 ±0.77 


16.09 ± 0.65 


16.17 ±0.68 


Rdist (-R©) 


16.42 ± 0.84 


16.42 ±0.84 


16.42 ±0.84 


16.42 ± 0.84 


16.50 ±0.84 


logg (cm s~ 2 ) 


3.394 ±0.016 


3.427 ± 0.008 


3.306 ±0.018 


3.302 ±0.012 


3.303 ± 0.018 


V^ot sini (km s _1 ) 


106.52 ± 6.39 


92.57 ±5.59 


79.91 ±4.79 


92.99 ±5.89 


92.99 ±4.62 


Xij 


26.11 


24.74 


21.14 


17.76 




xl 


46.71 


43.57 


24.96 


19.33 




f 


34.81 


32.49 


24.58 


23.99 




Xnv 


554.13 


551.57 


531.83 


536.32 




Xtot 


668.03 


652.54 


614.769 


597.67 




d.o.f. 


584 


583 


582 


581 





Note. — The assumed stellar temperature is T e g = 31,000 K. R op t is the stellar radius derived from 
the models. -R^ist is the stellar radius derived from the measured parallax and assumed temperature. 
V ro tsini = 95 ± 6 km s _1 is the projected rotational velocity of the O-star derived from the models. 
Model A assumes a circular orbit and synchronous rotation. Model B assumes a circular orbit and 
nonsynchronous rotation. Model C assumes an eccentric orbit and pseudo-synchronous rotation. Model 
D assumes an eccentric orbit and nonsynchronous rotation. The adopted parameters are the weighted 
averages of the values for model D derived for each temperature in the range of 30, 000 < T c ff < 32, 000 
K. 
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Table 2. Cygnus X-l Final Parameters 



parameter 


value 


% (deg) 


27.06 ± 0.76 


n 


1.400 ± 0.084 


e 


0.018 ±0.003 


u (deg) 


307.6 ±5.3 


M opt (Mq) 


19.16 ± 1.90 


Ropt (Rq) 


16.17 ±0.68 


log g (cgs) 


3.303 ±0.018 


M(M ) 


14.81 ±0.98 
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Temperature (kK) 



Fig. 1. — The derived O-star radius (top) and luminosity (bottom) as a function of the 
assumed effective temperature. 
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Fig. 2. — Top: The optical light curves and best-fitting models assuming an eccentric orbit 
with e = 0.018 (model D, left panels) and a circular orbit (model A, right panels). Note how 
much better the unequal maxima of the light curves are accommodated by the model that 
includes eccentricity. Bottom: The radial velocity measurements binned into 50 bins (filled 
circles) and the best-fitting model for the eccentric orbit (left) and the circular orbit (right). 
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Fig. 3. — A schematic diagram of Cygnus X-l, shown as it would appear on the sky plane. 
The offsets are in milliarcseconds (mas), assuming a distance of 1.86 kpc. The orbital phase is 
= 0.5, which corresponds to the superior conjunction of the O-star. The orbit of the black 
hole is indicated by the ellipse, where the major and minor axes have been drawn in for clarity 
(solid line and dashed line, respectivel y). The direction of the orbital motion is clockwise, as 
determined by the radio observations (iReid et al.l 1201 lh . The color map represents the local 
effective temperature. The star is much cool er near the inner Lagrang ian point because of 
the well-known effect of "gravity darkening" (jOrosz &: Hauschildtl 120001 ) . The temperatures 
referred to in the main text, Figure 1 and Table 1 refer to intensity- weighted average values. 
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A. ELC Fitting Details 

We computed models for 14 values of T eS between 28,000 and 34,000 K. We fit for each 
temperature separately since the derived radius from the parallax distance depends on the 
temperature, and the derived radius is used as an external constraint. For each model at 
each temperature, ELC's genetic code was run twice using different initial populations and 
the Monte Carlo Markov Chain (MCMC) code was run once. The best solutions were then 
refined using a simple grid search. We computed uncertainties on the fitted parameters and 
on the derived parameters (e.g. the black hol e mass M, the gravity of the O-star logg, etc.) 



using the procedure discussed in lOrosz et al.1 (120021 ) . The results for all 14 temperatures are 



shown in Table Al and Figs. Al and A2. As discussed in the main text, the improvement in 
the x 2 values as one goes from Model A to Model D is evident. For Model A, we consistently 
find -R op t < -Rdist- Furthermore, the rotational velocity derived from the model is consistently 
larger than the observed value. By allowing nonsynchronous rotation for the circular orbit 
(as in Model B), the model-derived stellar radius agrees with the radius computed from the 
distance, and the model-derived rotational velocity agrees with the measured one. Generally 
speaking, the star rotates slower than its synchronous value. However, an inspection of the 
light curves (see Figure 2 in the main text) shows that the maximum near phase 0.25 is 
slightly higher than the maximum near 0.75. Because an ellipsoidal model predicts maxima 
of equal intensity, the fit to the data is not optimal. We accommodate the unequal maxima 
by adding eccentricity to the synchronous model (as in Model C). However, in this case, 
the model stellar radius R opt is consistently larger than the radius derived from the distance 
i?dist and the model rotational velocity is smaller than the observed value. By allowing 
nonsynchronous rotation in the eccentric orbit (Model D), R opt agrees with -Rdist and the 
model rotational velocity agrees with the measured value. 



Table Al. Cygnus X-l Model Parameters 



T e ff 


model 


i 


n 


e 




M opt 


M 


Ropt a 




log g 


Vrot sin i c 


Xn 


Xr 






V 2 


(K) 




(deg) 






(deg) 


(M Q ) 


(Mq) 


(Rq) 


(Rq) 


(cgs) 


(km s~ 1 ) 












28000 


A d 


48.96 


1.000 






26.61 


9.01 


16.88 


19.26 


3.408 


115.08 


31.13 


59.68 


41.26 


556.00 


705.14 






±2.58 








±2.81 


±0.98 


±0.82 


±0.98 


±0.009 


±4.05 












28000 


B° 


65.50 


0.573 






40.56 


9.31 


19.85 


19.26 


3.450 


93.49 


27.53 


49.29 


35.33 


551.63 


664.20 






±1.59 


±0.039 






±4.14 


±0.58 


±0.85 


±0.98 


±0.008 


±5.97 












28000 


C f 


26.79 


1.000 


0.022 


296.4 


31.59 


19.58 


20.33 


19.26 


3.321 


86.79 


18.85 


23.13 


28.37 


533.37 


606.77 






±0.34 




±0.002 


±3.6 


±3.93 


±1.26 


±0.98 


±0.98 


±0.015 


±4.87 












28000 


Ds 


27.17 


1.187 


0.021 


298.1 


28.76 


18.31 


19.30 


19.26 


3.325 


94.80 


18.94 


22.83 


27.43 


534.25 


603.45 






±0.33 


±0.116 


±0.002 


±3.9 


±3.74 


±1.37 


±1.24 


±0.98 


±0.012 


±5.28 












28500 


A 


48.96 


1.000 






25.45 


8.77 


16.55 


18.53 


3.406 


112.83 


30.69 


58.09 


40.45 


555.65 


698.08 






±2.55 








±2.25 


±0.74 


±0.59 


±0.95 


±0.009 


±4.20 












28500 


B 


66.79 


0.591 






36.48 


8.62 


18.87 


18.53 


3.448 


92.64 


27.13 


48.42 


35.25 


551.91 


662.98 






±1.69 


±0.041 






±3.07 


±0.54 


±0.73 


±0.95 


±0.008 


±5.64 












28500 


C 


26.66 


1.000 


0.022 


296.4 


29.36 


18.89 


19.77 


18.53 


3.314 


84.02 


18.81 


23.04 


27.89 


533.17 


607.90 






±0.25 




±0.002 


±4.1 


±2.98 


±1.11 


±0.73 


±0.95 


±0.010 


±3.13 












28500 


D 


27.14 


1.263 


0.020 


299.6 


25.18 


17.04 


18.19 


18.53 


3.319 


94.96 


18.75 


22.29 


26.76 


534.55 


602.38 






±0.28 


±0.112 


±0.002 


±4.6 


±3.40 


±1.27 


±1.06 


±0.95 


±0.010 


±5.28 












29000 


A 


48.91 


1.000 






23.88 


8.45 


16.12 


17.88 


3.401 


109.80 


29.77 


55.81 


39.34 


555.14 


689.30 






±2.45 








±1.83 


±0.78 


±0.64 


±0.99 


±0.010 


±3.97 












29000 


B 


67.21 


0.626 






32.06 


7.93 


17.82 


17.88 


3.442 


92.97 


26.76 


47.33 


34.66 


551.77 


660.64 






±2.30 


±0.053 






±5.35 


±0.94 


±1.31 


±0.99 


±0.009 


±5.83 












29000 


C 


26.59 


1.000 


0.022 


297.0 


28.37 


18.60 


19.52 


17.88 


3.310 


82.70 


18.64 


22.72 


27.48 


533.05 


608.72 






±0.31 




±0.002 


±5.0 


±2.40 


±0.81 


±0.63 


±0.99 


±0.010 


±2.74 












29000 


D 


27.13 


1.330 


0.020 


301.5 


22.28 


15.95 


17.24 


17.88 


3.313 


94.70 


18.30 


21.18 


25.67 


535.31 


601.21 






±0.35 


±0.110 


±0.002 


±5.2 


±2.84 


±1.15 


±0.96 


±0.99 


±0.010 


±4.97 












29500 


A 


48.85 


1.000 






23.48 


8.37 


16.02 


17.34 


3.399 


109.01 


28.89 


53.39 


38.06 


555.08 


683.10 






±2.45 








±2.19 


±0.60 


±0.49 


±0.89 


±0.012 


±3.86 












29500 


B 


67.66 


0.664 






28.02 


7.27 


16.79 


17.34 


3.435 


93.19 


26.61 


46.69 


34.43 


551.80 


660.02 






±2.88 


±0.040 






±3.12 


±0.53 


±0.83 


±0.89 


±0.006 


±5.74 












29500 


C 


29.00 


1.000 


0.028 


300.7 


25.80 


15.72 


18.45 


17.34 


3.317 


85.86 


22.03 


27.08 


26.37 


532.01 


611.39 






±0.40 




±0.002 


±4.1 


±0.62 


±0.32 


±0.17 


±0.89 


±0.005 


±1.08 












29500 


D 


27.13 


1.347 


0.019 


302.4 


21.32 


15.57 


16.92 


17.34 


3.310 


94.20 


18.18 


20.80 


25.37 


535.26 


599.84 






±0.33 


±0.063 


±0.002 


±9.99 


±1.73 


±0.79 


±0.63 


±0.89 


±0.011 


±4.70 













Table Al — Continued 





model 


i 


n 


e 




M opt 


M 


-Ropt a 


R b 


logfl 


Vrot sin i c 


xl 




X 2 y 




Xtot 


(K) 




(deg) 






(deg) 


(M@) 


(M ) 


(Re) 


(Re) 


(cgs) 


(km s- 1 ) 












29750 


A 


48.82 


1.000 






22.03 


8.06 


15.61 


17.11 


3.394 


106.17 


28.71 


52.94 


37.68 


554.57 


680.32 






±5.14 








±2.25 


±0.97 


±0.73 


±0.88 


±0.019 


±7.00 












29750 


B 


67.78 


0.675 






27.08 


7.11 


16.53 


17.11 


3.434 


93.29 


26.06 


45.74 


34.43 


551.80 


658.56 






±2.40 


±0.055 






±3.35 


±0.57 


±0.83 


±0.88 


±0.011 


±5.63 












29750 


C 


28.93 


1.000 


0.027 


301.5 


25.45 


15.65 


18.38 


17.11 


3.315 


85.26 


26.24 


26.82 


25.83 


531.92 


611.22 






±2.83 




±0.007 


±8.3 


±2.72 


±2.91 


±1.09 


±0.88 


±0.025 


±9.45 












29750 


D 


27.09 


1.358 


0.019 


303.4 


20.66 


15.34 


16.71 


17.11 


3.307 


93.56 


18.00 


20.13 


24.75 


535.50 


598.65 






±0.32 


±0.080 


±0.002 


±5.1 


±2.76 


±1.08 


±0.92 


±0.88 


±0.011 


±4.88 












30000 


A 


48.82 


1.000 






22.61 


8.19 


15.78 


16.91 


3.396 


107.36 


28.06 


51.49 


37.01 


554.77 


677.29 






±3.07 








±1.81 


±0.78 


±0.64 


±0.86 


±0.012 


±5.96 












30000 


B 


67.86 


0.683 






26.16 


6.96 


16.30 


16.91 


3.431 


93.17 


25.81 


45.18 


33.70 


551.81 


657.10 






±2.87 


±0.045 






±3.04 


±0.63 


±0.83 


±0.86 


±0.008 


±5.39 












30000 


C 


28.90 


1.000 


0.027 


302.4 


25.57 


15.71 


18.42 


16.91 


3.315 


85.33 


18.33 


26.29 


25.33 


531.90 


610.91 






±0.32 




±0.006 


±4.2 


±2.34 


±1.19 


±0.53 


±0.86 


±0.027 


±9.52 












30000 


D 


27.11 


1.376 


0.019 


304.3 


20.21 


15.14 


16.53 


16.91 


3.307 


93.86 


17.98 


20.19 


24.89 


535.63 


598.92 






±0.37 


±0.094 


±0.002 


±4.8 


±1.80 


±0.78 


±0.64 


±0.86 


±0.011 


±5.17 












30500 


A 


48.78 


1.000 






21.94 


8.04 


15.62 


16.61 


3.392 


106.19 


27.18 


49.42 


35.84 


554.54 


671.84 






±2.77 








±2.17 


±0.85 


±0.75 


±0.85 


±0.015 


±5.44 












30500 


B 


67.80 


0.681 






26.18 


6.96 


16.35 


16.61 


3.429 


93.19 


25.29 


44.34 


33.15 


551.76 


654.73 






±0.67 


±0.040 






±1.96 


±0.34 


±0.53 


±0.85 


±0.007 


±6.17 












30500 


C 


28.78 


1.000 


0.026 


303.5 


25.44 


15.75 


18.41 


16.61 


3.313 


84.79 


18.20 


25.31 


24.71 


531.95 


610.76 






±0.36 




±0.002 


±4.9 


±1.98 


±0.32 


±0.09 


±0.85 


±0.005 


±1.03 












30500 


D 


27.07 


1.392 


0.018 


306.3 


19.54 


14.90 


16.30 


16.61 


3.304 


93.46 


17.84 


19.71 


24.41 


536.01 


598.17 






±0.29 


±0.071 


±0.003 


±5.3 


±1.96 


±0.71 


±0.71 


±0.85 


±0.013 


±4.15 












31000 


A 


51.45 


1.000 






20.53 


7.37 


15.07 


16.42 


3.394 


106.52 


26.11 


46.71 


34.81 


554.13 


668.03 






±7.46 








±2.00 


±1.19 


±1.26 


±0.84 


±0.016 


±6.39 












31000 


B 


67.74 


0.674 






26.27 


6.98 


16.41 


16.42 


3.427 


92.57 


24.74 


43.57 


32.49 


551.57 


652.54 






±2.83 


±0.043 






±2.41 


±0.39 


±0.72 


±0.84 


±0.008 


±5.59 












31000 


C 


28.50 


1.000 


0.026 


303.5 


25.17 


15.83 


18.46 


16.42 


3.306 


79.91 


21.14 


24.96 


24.58 


531.83 


614.77 






±2.21 




±0.003 


±5.1 


±2.54 


±1.80 


±0.77 


±0.84 


±0.018 


±4.79 












31000 


D 


27.03 


1.404 


0.018 


308.1 


18.97 


14.69 


16.09 


16.42 


3.302 


92.99 


17.76 


19.33 


23.99 


536.32 


597.67 






±0.41 


±0.099 


±0.002 


±5.5 


±2.15 


±0.70 


±0.65 


±0.84 


±0.012 


±5.89 













Table Al — Continued 



T c tt 


model 


i 


n 


e 




M opt 


M 


d a 
-n-opt 


-Rdist b 


logs 


Vrot sin i c 


xl 




X 2 y 


xh 


Xtot 


(K) 




(deg) 






(deg) 


(M @ ) 


(M Q ) 


(Re) 


(Re) 


(cgs) 


(km s- 1 ) 












31500 


A 


46.47 


1.000 






20.94 


8.19 


15.53 


16.31 


3.376 


101.78 


26.11 


47.59 


34.20 


554.08 


664.12 






±2.89 








±2.29 


±0.77 


±0.71 


±0.83 


±0.016 


±5.21 












31500 


B 


67.64 


0.673 






26.31 


6.99 


16.49 


16.31 


3.424 


92.74 


24.34 


43.24 


31.98 


551.46 


651.22 






±2.95 


±0.045 






±3.10 


±0.39 


±0.80 


±0.83 


±0.011 


±6.18 












31500 


C 


28.67 


1.000 


0.025 


306.1 


24.75 


15.58 


18.26 


16.31 


3.308 


83.55 


20.90 


24.24 


23.62 


532.01 


609.88 






±0.40 




±0.004 


±6.3 


±2.55 


±0.98 


±0.75 


±0.83 


±0.014 


±4.16 












31500 


D 


26.96 


1.412 


0.017 


310.0 


18.56 


14.57 


15.95 


16.31 


3.301 


92.43 


17.69 


19.02 


23.59 


536.64 


597.21 






±1.30 


±0.078 


±0.005 


±5.5 


±1.83 


±1.27 


±0.65 


±0.83 


±0.012 


±3.80 












32000 


A 


46.38 


1.000 






20.73 


8.16 


15.51 


16.26 


3.373 


101.49 


25.76 


47.03 


33.69 


553.80 


662.26 






±4.79 








±2.02 


±0.93 


±0.70 


±0.83 


±0.017 


±6.02 












32000 


B 


67.96 


0.703 






23.60 


6.53 


15.74 


16.26 


3.417 


92.72 


24.05 


42.77 


31.82 


551.29 


650.46 






±3.06 


±0.047 






±2.56 


±0.34 


±0.93 


±0.83 


±0.007 


±6.11 












32000 


C 


28.60 


1.000 


0.024 


306.0 


24.78 


15.64 


18.29 


16.26 


3.307 


83.42 


20.68 


23.77 


23.34 


532.13 


609.55 






±1.60 




±0.002 


±6.8 


±0.52 


±1.46 


±0.45 


±0.83 


±0.016 


±3.01 












32000 


D 


26.97 


1.411 


0.017 


310.0 


18.56 


14.56 


15.95 


16.26 


3.301 


92.43 


17.70 


18.91 


23.48 


536.62 


597.03 






±1.44 


±0.076 


±0.005 


±5.3 


±1.74 


±1.43 


±0.75 


±0.83 


±0.011 


±4.11 












32500 


A 


46.32 


1.000 






20.58 


8.13 


15.50 


16.26 


3.371 


101.30 


25.41 


46.45 


33.29 


553.55 


660.64 






±5.54 








±2.26 


±1.32 


±0.92 


±0.84 


±0.019 


±7.42 












32500 


B 


67.91 


0.705 






23.44 


6.50 


15.74 


16.26 


3.414 


92.91 


23.66 


42.52 


31.65 


551.14 


649.47 






±7.48 


±0.058 






±2.67 


±0.46 


±0.75 


±0.84 


±0.013 


±5.32 












32500 


C 


28.56 


1.000 


0.024 


306.0 


24.97 


15.74 


18.36 


16.26 


3.307 


83.58 


20.58 


23.60 


23.34 


532.18 


609.61 






±0.43 




±0.002 


±5.4 


±0.27 


±0.26 


±0.21 


±0.84 


±0.005 


±1.36 












32500 


D 


26.99 


1.412 


0.017 


310.7 


18.57 


14.55 


15.96 


16.26 


3.301 


92.57 


17.81 


19.01 


23.21 


536.58 


596.90 






±0.93 


±0.091 


±0.003 


±6.1 


±2.38 


±0.73 


±0.76 


±0.84 


±0.011 


±5.01 












33000 


A 


46.30 


1.000 






20.50 


8.12 


15.50 


16.29 


3.369 


101.27 


25.13 


46.01 


32.99 


553.36 


659.47 






±4.84 








±2.42 


±1.44 


±0.92 


±0.84 


±0.020 


±5.90 












33000 


B 


67.76 


0.693 






24.18 


6.63 


15.99 


16.29 


3.413 


92.76 


23.57 


42.71 


31.39 


551.05 


648.98 






±5.30 


±0.049 






±2.07 


±0.38 


±0.61 


±0.84 


±0.008 


±5.89 












33000 


C 


28.50 


1.000 


0.024 


306.9 


25.06 


15.82 


18.40 


16.29 


3.307 


83.59 


20.66 


23.86 


22.96 


532.20 


609.65 






±0.42 




±0.002 


±5.8 


±0.26 


±0.31 


±0.07 


±0.84 


±0.005 


±1.15 












33000 


D 


27.89 


1.386 


0.020 


312.0 


19.87 


14.42 


16.34 


16.29 


3.310 


95.95 


18.56 


20.65 


22.32 


535.30 


596.86 






±1.15 


±0.077 


±0.005 


±5.7 


±2.27 


±1.20 


±0.90 


±0.84 


±0.016 


±6.65 













Table Al — Continued 





model 


i 


n 


e 




M opt 


M 


D a 

JLopt 


^dist b 


logfl 


Vrot sini c 


xl 


xl 


xl 


Xrv 


Xtot 


(K) 




(deg) 






(deg) 


(M@) 


(M e ) 


(Re) 


(Re) 


(cgs) 


(km s _1 ) 












33500 


A 


44.04 


1.000 






20.66 


8.57 


15.72 


16.32 


3.360 


98.75 


24.97 


46.15 


32.96 


553.02 


658.01 






±2.85 








±2.06 


±1.10 


±0.63 


±0.84 


±0.013 


±5.62 












33500 


13 


60.52 


0.730 






24.15 


7.13 


16.21 


16.32 


3.401 


93.14 


23.45 


42.50 


31.29 


551.31 


648.66 






±7.31 


±0.077 






±2.54 


±0.76 


±0.76 


±0.84 


±0.014 


±5.52 












33500 


C 


28.49 


1.000 


0.024 


305.9 


25.07 


15.83 


18.42 


16.32 


3.306 


83.62 


20.61 


24.09 


23.06 


532.27 


609.97 






±2.16 




±0.004 


±6.1 


±1.06 


±1.09 


±0.22 


±0.84 


±0.019 


±5.56 












33500 


D 


27.81 


1.398 


0.019 


312.0 


19.70 


14.42 


16.28 


16.32 


3.309 


96.13 


18.56 


20.64 


22.30 


535.26 


596.80 






±0.87 


±0.079 


±0.003 


±5.5 


±2.35 


±1.29 


±0.92 


±0.84 


±0.014 


±5.39 












34000 


A 


44.04 


1.000 






20.63 


8.56 


15.73 


16.34 


3.359 


98.81 


24.73 


45.69 


32.69 


552.78 


656.84 






±2.80 








±1.64 


±0.95 


±0.66 


±0.84 


±0.013 


±5.24 












34000 


B 


65.32 


0.703 






24.04 


6.75 


16.08 


16.34 


3.406 


92.80 


23.30 


42.94 


31.24 


550.83 


648.55 






±7.42 


±0.085 






±2.36 


±0.50 


±0.64 


±0.84 


±0.016 


±5.32 












34000 


C 


28.38 


1.000 


0.024 


306.0 


25.06 


15.90 


18.43 


16.34 


3.306 


83.37 


22.29 


24.31 


23.01 


532.23 


610.29 






±2.32 




±0.005 


±5.1 


±0.79 


±1.41 


±0.55 


±0.84 


±0.017 


±5.88 












34000 


D 


27.83 


1.405 


0.019 


312.0 


19.71 


14.41 


16.27 


16.34 


3.310 


96.62 


18.62 


20.78 


22.38 


535.19 


597.05 






±0.97 


±0.073 


±0.003 


±5.8 


±1.94 


±0.61 


±0.60 


±0.84 


±0.013 


±5.29 













a The stellar radius derived from the models. 

b The stellar radius derived from the measured parallax and the assumed temperature. 

c The projected rotational velocity of the O-star derived from the models. The observed value is 95 ± 6 km s -1 . 

d Model A assumes a circular orbit and synchronous rotation. 

c Modcl B assumes a circular orbit and nonsynchronous rotation. 

f Model C assumes an eccentric orbit and pseudo-synchronous rotation 

g Model D assumes an eccentric orbit and nonsynchronous rotation 




28 29 30 31 32 33 34 



T eff ( kK ) 

Fig. Al. — Results of the dynamical modeling assuming a circular orbit. Various quantities 
are plotted as a function of temperature for the model assuming synchronous rotation (open 
circles) and nonsynchronous rotation (filled circles). From top to bottom the plots show 
(a) the total x 2 ; (b) the difference between R dlst and R 2 in solar radii; (c) the difference 
between the observed projected rotational velocity and the value derived from the models 
in km s _1 , where the hatched region denotes the la uncertainty on the measured value; (d) 
the inclination in degrees; (e) the mass of the O-star in solar masses; and (f) the mass of the 
black hole in solar masses. In general, for the synchronous models the O-star radius derived 
from the models is smaller than the radius derived from the distance, whereas the derived 
projected rotational velocity is larger than observed. 
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Fig. A2. — Same as Figure IA14 but for the models that assume an eccentric orbit. In this 
case, the 0-star radius derived from the models is larger than the radius derived from the 
distance, and the derived rotational velocity is smaller than observed. 
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